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Abstract 

We have studied kinetics of random sequential adsorption of mix- 
tures on a square lattice using Monte Carlo method. Mixtures of linear 
short segments and long segments were deposited with the probabil- 
ity p and 1 — p, respectively. For fixed lengths of each segment in the 
mixture, the jamming limits decrease when p increases. The jamming 
limits of mixtures always are greater than those of the pure short- 
or long-segment deposition. For fixed p and fixed length of the short 
segments, the jamming limits have a maximum when the length of 
the long segment increases. We conjectured a kinetic equation for the 
jamming coverage based on the data fitting. 
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Random sequential adsorption(RSA) is a model of irreversible deposi- 
tion of fixed-shape objects[l-3]. The objects are deposited randomly and 
irreversibly on a substrate. When a depositing object overlaps a deposited 
object, the depositing object is removed from the system and another de- 
positing is attempted sequentially. When we allow a formation of a single 
layer, the substrate does not reach the full coverage at the long time limit. 
The jamming coverage 6(t) is defined by the total number of covered sites 
divided by the lattice size. The coverage converges to a particular value, 
jamming limit 0(oo), in the long time limit. RSA is realized in experimen- 
tal studies of protein and colloid particle adhesion on surfaces under the 
conditions of negligibly slow surface relaxation[4-7]. 

The theoretical studies including rate equations [8- 10], series expansionfll], 
and Monte Carlo method[12-18] were reported for both continuous and dis- 
crete models. Exact solutions are also available for one-dimensional mod- 
els[13, 14,19-21]. In lattice deposition models the jamming coverage is asymp- 
totically exponential 

0(f) = 0(oo) - Aexp(-Bt) (1) 

where A and B are parameters which depend on the dimensionality of the 
substrate and on the shape of the depositing objects. There are several stud- 
ies of kinetics of mixture depositions [22-28]. Two different kinds of objects 
are deposited on the substrate with a different adsorption probability of each 
object. Exact solutions were reported for RSA of mixtures of monomer and 
linear /c-mers in a one-dimensional lattice [28]. They showed that the ad- 
dition of pointlike particles modifies in a nonuniversal way the form of the 
long-time convergence law of the approach to the jamming coverage. RSA of 
arbitrary mixtures of line segments of two different lengths were solved an- 
alytically on the one-dimensional lattice[29-30]. The deposition on a square 
lattice of a mixture of line segments of length k\ and k 2 (end-on model), 
chosen with equal probability, have been studied by Svrakic and Henkel us- 
ing Monte Carlo method[26]. They found two inequality relations for the 
jamming limits 

0(fci,fc 2 ,oo) > 0(fci,oo) > 9(k 2 ,oo) for k 2 > h (2) 

and 

0(&4, k' 2 , oo) > 0(&4, k 2 , oo) for k' 2 >k\ (3) 
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where 9(k\, k 2 , oo) is the jamming limit of the jamming coverage 9(ki, k 2 ,t) 
of a mixture and 9(k, oo) is the jamming limit of single k-mer. They also 
proposed the time dependence of the jamming coverage 

9(h, k 2 , t) = 9(h, k 2 , oo) - A(k 1 , k 2 ) exp(-B(k u k 2 )t) (4) 

with B(ki, k 2 ) ~ 1.0 for k\ ^ k 2 . 

In the present work we have studied kinetics of RSA of mixture having a 
general deposition probability of each segment, that is, the deposition prob- 
ability p for short length of linear /comers and 1 — p for long length of linear 
/c 2 -mers. Monte Carlo simulations have been performed on a square lattice 
L x L. We used lattices of size L = 512. We randomly select a lattice site 
and try to deposit a line segment of length k\ with the probability p or length 
k 2 with the probability 1 — p. We always take k\ < k 2 . If a chosen site is 
occupied by a deposited object, the attempt fails, the time is increased by 
one unit, and a new site and a new line segment are selected. If a chosen 
site is empty, then we randomly select a direction of the four possible ori- 
entations. If all k± — 1 (or k 2 — 1) neighbor sites are empty, the selected 
sites are occupied. If any neighbor site is occupied by previously deposited 
object in chosen direction, the attempt fails and the time is increased by 
one unit and we try a new site. This model is standard model of RSA [25]. 
In end-on model, all possible directions of a selected site are consumed by 
occuping selected line segments. After a long time, the system is close to 
the jamming limit. At late-time stage we check all empty sites. If an empty 
site has at least one possible direction which is available to accummulate a 
short line segment fci-mer, then empty site is marked as an accessible site. 
If an empty site has no possible direction to occupy fci-mer, that site is an 
inaccessible site. We further try to occupy line segments at accessible sites. 
If there are no accessible sites, the system is in the jamming limit. We always 
use periodic boundary conditions. One Monte Carlo time step is defined by 
the total number of attempts to select a site divided by the total number of 
lattice sites. The data are averaged over 50 independent runs for each choice 
of mixtures and each deposition probability. 

We present the typical jamming configuration of mixtures (k\ — 4, k 2 — 
16) in Fig. 1 (a) and {k\ — 4, k 2 — 32) in Fig. 1 (b) for the short A^-mer 
deposition probability p = 0.5. The lattice size is 128 x 128 square lattice 
with periodic boundary conditions. The configurational structure of jamming 
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limit has a lot of local structures which have parallel short and long line 
segments. In Fig. 2 we plot \n[6(ki, k 2 , oo)—6(ki, k 2 , t)] versus time for k± — 2 
and p — 0.1. At long time limit the lines are linear and parallel. In such a 
standard model the jamming coverages follow the exponential behaviour like 
eq. (4). We have calculated the slopes of the parallel lines for various length 
of (ki,k 2 ) pairs and p. We concluded B(ki, k 2 ,p) = p/2 from linear slopes. 
In Fig. 2 we can observe the independence of B(ki,k 2 ,p) on the length 
of each line segment. All lines are parallel and the slopes depend only on 
the deposition probability. These results are different from the observation 
of end-on model B(k\,k 2 ,p = 1/2) = 1 [26] and single adsorption model 
B(k) = 2[31]. In end-on model we observed B(ki, k 2 ,p) = 2p[32]. 

In Fig. 3 we show the jamming limits versus the short fci-mer deposi- 
tion probability p for (ki = 2, k 2 ) mixtures. The jamming limits decrease 
monotonically when p increases. The jamming limit approaches to a value 
of 8(k = 2, oo) = 0.9068 at p — >• 1. At p — > the jamming limit shows a 
singularity. For example, the jamming limit drops to a value of pure single 
segment deposition 6{k\ — 2, k 2 — 3, oo) = 0.8467 at p — 0. This singu- 
larity was already observed in one-dimensional adsorption of mixtures [30]. 
For fixed length of k\ and k 2 we observed that the jamming limit satisfied 
the same inequality relation as eq. (2). The mixture depositions lead to 
more efficient jamming coverage than depositions of the single line segment. 
The jamming coverages are controlled by the long segment of mixtures up 
to cross-over times t x . At t > t x . The local empty spaces are shorter than 
the length of the long segments. Therefore, the short segments are further 
adsorbed on the substrate. 

In Fig. 4 we plot the jamming limit 9(k\ = 2, k 2 , oo;p) versus the length 
of a long segment for fixed length of short segment k\ — 2 and various p- 
values. For a fixed value k 2 , the jamming limits satisfy 9(ki,k 2 ,oo;p 2 ) > 
9(ki, k 2 , oo;pi) when p 2 > p 1 . For a fixed value p the jamming limits reach 
a maximum value at a particular value k%- A combination (A^,/^) of mix- 
tures at a fixed value of p gives the most efficient jamming coverage at long 
time limit. We check the finite size effects for p = 0.5 using 256 x 256 and 
1024 x 1024 lattices. We observed that the effects of finite size are negligible. 
We concluded that these maximum behaviours are intrinsic in mixture depo- 
sition on such a standard model. We also observed the maximum behaviour 
for end-on model at a certain range of the probability p. (These results will 
be published elsewhere [32].) When the length of the long segments increases 
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for fixed length of ki, the adsorbed long segments induce that the newly 
adsorbing long segments locally have the same direction as adsorbed one as 
shown in Fig. 1 (b). Among the absorbed parallel long segments, the short 
A^-mers further adsorbed to the direction of the long /c 2 -mers. These adsorp- 
tions induced local one-dimensional characteristics of the short segments and 
the jamming limits decrease for fc 2 > ^ 2 - 

We have performed the data fitting to check the parameter dependence 
of the amplitude A(k±, k 2 ,p). We concluded that the amplitude has the 
functional form as A(ki,k2,p) = C (l//ci) 2 exp[Ci(l — p) + C 2 (l —p)/k?\. 
The amplitudes are well fitted for the choice of the constant C Q = 0.8 ±0.1, 
C± = — 0.5±0.2 and C 2 = — 1.3±0.1 for all k\ < fc 2 mixture combinations [32]. 

In summary we observed that the mixture depositions cover the substrate 
more efficiently than the deposition of single-length segments. We observed 
that the jamming limits show a maximum behaviour for fixed k\ and p when 
the lengths of fc 2 vary. The jamming coverage shows exponential behaviour. 
We proposed a functional form of the amplitude of coverage from the data 
fitting. 

This work has been supported by Inha University and by the Basic Science 
Institute Program, Ministry of Education, Project No. BSRI-97-2430. I 
wish to thank Professor Vladimir Privman for his critical reading of this 
manuscript. 
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Figure Captions 



Figure 1: Typical configurations of mixture deposition on 128 x 128 square 
lattices for (a) ki = A,k 2 = 16 and (b) ki = 4, k 2 = 32. 

Figure 2: The plot of \n[9(ki, k 2 , oo) — 9(k±, k 2 , t)] versus time for k\ = 2, 
p — 0.1 and k 2 —4, 8, 16, 32, 64 from bottom to top. 

Figure 3: The jamming limits versus the short segment deposition probability 
p for k 1 = 2 and k 2 = 3(o), 4(+), 8(*), 16(A), 32(D), and 64(x). 

Figure 4: The jamming limits 6{k\ = 2, k 2 , oo; p) versus the length of the long 
segment for a fixed length of the short segment hi — 2 and p = 0.1(A), 
0.3(x), 0.5(D), 0.7(+), and 0.9(o). 
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